Cell-type-specific optogenetic fMRI on basal forebrain reveals functional network basis of behavioral preference.

The basal forebrain (BF) is a complex structure that plays key roles in regulating various brain functions. However, it remains unclear how cholinergic and non-cholinergic BF neurons modulate large-scale functional networks and their relevance in intrinsic and extrinsic behaviors. With an optimized awake mouse optogenetic fMRI approach, we revealed that optogenetic stimulation of four BF neuron types evoked distinct cell-type-specific whole-brain BOLD activations, which could be attributed to BF-originated low-dimensional structural networks. Additionally, optogenetic activation of VGLUT2, ChAT, and PV neurons in the BF modulated the preference for locomotion, exploration, and grooming, respectively. Furthermore, we uncovered the functional network basis of the above BF-modulated behavioral preference through a decoding model linking the BF-modulated BOLD activation, low-dimensional structural networks, and behavioral preference. To summarize, we decoded the functional network basis of differential behavioral preferences with cell-type-specific optogenetic fMRI on the BF and provided an avenue for investigating mouse behaviors from a whole-brain view.

Differential Effect of Global Signal Regression Between Awake and Anesthetized Conditions in Mice.

Introduction: In resting-state functional magnetic resonance imaging (rs-fMRI) studies, global signal regression (GSR) is a controversial preprocessing strategy. It effectively eliminates global noise driven by motion and respiration but also can introduce artifacts and remove functionally relevant metabolic information. Most preclinical rs-fMRI studies are performed in anesthetized animals, and anesthesia will alter both metabolic and neuronal activity. Methods: In this study, we explored the effect of GSR on rs-fMRI data collected under anesthetized and awake state in mice (n = 12). We measured global signal amplitude, and also functional connectivity (FC), functional connectivity density (FCD) maps, and brain modularity, all commonly used data-driven analysis methods to quantify connectivity patterns. Results: We found that global signal amplitude was similar between the awake and anesthetized states. However, GSR had a different impact on connectivity networks and brain modularity changes between states. We demonstrated that GSR had a more prominent impact on the anesthetized state, with a greater decrease in functional connectivity and increased brain modularity. We classified mice using the change in amplitude of brain modularity coefficient (ΔQ) before and after GSR processing. The results revealed that, when compared with the largest ΔQ group, the smallest ΔQ group had increased FCD in the cortex region in both the awake and anesthetized states. This suggests differences in individual mice may affect how GSR differentially affects awake versus anesthetized functional connectivity. Discussion: This study suggests that, for rs-fMRI studies which compare different physiological states, researchers should use GSR processing with caution.

State-dependent memory retrieval: insights from neural dynamics and behavioral perspectives.

Memory retrieval is strikingly susceptible to external states (environment) and internal states (mood states and alcohol), yet we know little about the underlying mechanisms. We examined how internally generated states influence successful memory retrieval using the functional magnetic resonance imaging (fMRI) of laboratory mice during memory retrieval. Mice exhibited a strong tendency to perform memory retrieval correctly only in the reinstated mammillary body-inhibited state, in which mice were trained to discriminate auditory stimuli in go/no-go tasks. fMRI revealed that distinct auditory cues engaged differential brain regions, which were primed by internal state. Specifically, a cue associated with a reward activated the lateral amygdala, while a cue signaling no reward predominantly activated the postsubiculum. Modifying these internal states significantly altered the neural activity balance between these regions. Optogenetic inhibition of those regions in the precue period blocked the retrieval of type-specific memories. Our findings suggest that memory retrieval is under the control of two interrelated neural circuits underlying the neural basis of state-dependent memory retrieval.

An atlas of itch-associated neural dynamics in the mouse brain.

The itch-scratching cycle is mediated by neural dynamics in the brain. However, our understanding of the neural dynamics during this cycle remains limited. In this study, we examine the neural dynamics of 126 mouse brain areas by measuring the calcium signal using fiber photometry. We present numerous response patterns in the mouse brain during the itch-scratching cycle. Interestingly, we find that a group of brain areas exhibit activation only at the end of histamine-induced scratching behavior. Additionally, several brain areas exhibit transient activation at the onset of scratching induced by chloroquine. Both histamine- and chloroquine-induced itch evoke diverse response patterns across the mouse brain. In summary, our study provides a comprehensive dataset for the diverse activity pattern of mouse brain during the itch-scratching cycle, paving the way for further exploration into the neural mechanisms underlying the itch-scratching cycle.

Anesthesia-related brain microstructure modulations detected by diffusion magnetic resonance imaging.

Recent studies have shown significant changes to brain microstructure during sleep and anesthesia. In vivo optical microscopy and magnetic resonance imaging (MRI) studies have attributed these changes to anesthesia and sleep-related modulation of the brain's extracellular space (ECS). Isoflurane anesthesia is widely used in preclinical diffusion MRI (dMRI) and it is therefore important to investigate if the brain's microstructure is affected by anesthesia to an extent detectable with dMRI. Here, we employ diffusion kurtosis imaging (DKI) to assess brain microstructure in the awake and anesthetized mouse brain (n = 22). We find both mean diffusivity (MD) and mean kurtosis (MK) to be significantly decreased in the anesthetized mouse brain compared with the awake state (p < 0.001 for both). This effect is observed in both gray matter and white matter. To further investigate the time course of these changes we introduce a method for time-resolved fast DKI. With this, we show the time course of the microstructural alterations in mice (n = 5) as they transition between states in an awake-anesthesia-awake paradigm. We find that the decrease in MD and MK occurs rapidly after delivery of gas isoflurane anesthesia and that values normalize only slowly when the animals return to the awake state. Finally, time-resolved fast DKI is employed in an experimental mouse model of brain edema (n = 4), where cell swelling causes the ECS volume to decrease. Our results show that isoflurane affects DKI parameters and metrics of brain microstructure and point to isoflurane causing a reduction in the ECS volume. The demonstrated DKI methods are suitable for in-bore perturbation studies, for example, for investigating microstructural modulations related to sleep/wake-dependent functions of the glymphatic system. Importantly, our study shows an effect of isoflurane anesthesia on rodent brain microstructure that has broad relevance to preclinical dMRI.

Single-cell spatial transcriptome reveals cell-type organization in the macaque cortex.

Elucidating the cellular organization of the cerebral cortex is critical for understanding brain structure and function. Using large-scale single-nucleus RNA sequencing and spatial transcriptomic analysis of 143 macaque cortical regions, we obtained a comprehensive atlas of 264 transcriptome-defined cortical cell types and mapped their spatial distribution across the entire cortex. We characterized the cortical layer and region preferences of glutamatergic, GABAergic, and non-neuronal cell types, as well as regional differences in cell-type composition and neighborhood complexity. Notably, we discovered a relationship between the regional distribution of various cell types and the region's hierarchical level in the visual and somatosensory systems. Cross-species comparison of transcriptomic data from human, macaque, and mouse cortices further revealed primate-specific cell types that are enriched in layer 4, with their marker genes expressed in a region-dependent manner. Our data provide a cellular and molecular basis for understanding the evolution, development, aging, and pathogenesis of the primate brain.

Instantaneous antidepressant effect of lateral habenula deep brain stimulation in rats studied with functional MRI.

The available treatments for depression have substantial limitations, including low response rates and substantial lag time before a response is achieved. We applied deep brain stimulation (DBS) to the lateral habenula (LHb) of two rat models of depression (Wistar Kyoto rats and lipopolysaccharide-treated rats) and observed an immediate (within seconds to minutes) alleviation of depressive-like symptoms with a high-response rate. Simultaneous functional MRI (fMRI) conducted on the same sets of depressive rats used in behavioral tests revealed DBS-induced activation of multiple regions in afferent and efferent circuitry of the LHb. The activation levels of brain regions connected to the medial LHb (M-LHb) were correlated with the extent of behavioral improvements. Rats with more medial stimulation sites in the LHb exhibited greater antidepressant effects than those with more lateral stimulation sites. These results indicated that the antidromic activation of the limbic system and orthodromic activation of the monoaminergic systems connected to the M-LHb played a critical role in the rapid antidepressant effects of LHb-DBS. This study indicates that M-LHb-DBS might act as a valuable, rapid-acting antidepressant therapeutic strategy for treatment-resistant depression and demonstrates the potential of using fMRI activation of specific brain regions as biomarkers to predict and evaluate antidepressant efficacy.

Sleep fMRI with simultaneous electrophysiology at 9.4 T in male mice.

Sleep is ubiquitous and essential, but its mechanisms remain unclear. Studies in animals and humans have provided insights of sleep at vastly different spatiotemporal scales. However, challenges remain to integrate local and global information of sleep. Therefore, we developed sleep fMRI based on simultaneous electrophysiology at 9.4 T in male mice. Optimized un-anesthetized mouse fMRI setup allowed manifestation of NREM and REM sleep, and a large sleep fMRI dataset was collected and openly accessible. State dependent global patterns were revealed, and state transitions were found to be global, asymmetrical and sequential, which can be predicted up to 17.8 s using LSTM models. Importantly, sleep fMRI with hippocampal recording revealed potentiated sharp-wave ripple triggered global patterns during NREM than awake state, potentially attributable to co-occurrence of spindle events. To conclude, we established mouse sleep fMRI with simultaneous electrophysiology, and demonstrated its capability by revealing global dynamics of state transitions and neural events.

An integrated resource for functional and structural connectivity of the marmoset brain.

Comprehensive integration of structural and functional connectivity data is required to model brain functions accurately. While resources for studying the structural connectivity of non-human primate brains already exist, their integration with functional connectivity data has remained unavailable. Here we present a comprehensive resource that integrates the most extensive awake marmoset resting-state fMRI data available to date (39 marmoset monkeys, 710 runs, 12117 mins) with previously published cellular-level neuronal tracing data (52 marmoset monkeys, 143 injections) and multi-resolution diffusion MRI datasets. The combination of these data allowed us to (1) map the fine-detailed functional brain networks and cortical parcellations, (2) develop a deep-learning-based parcellation generator that preserves the topographical organization of functional connectivity and reflects individual variabilities, and (3) investigate the structural basis underlying functional connectivity by computational modeling. This resource will enable modeling structure-function relationships and facilitate future comparative and translational studies of primate brains.

Multimodal analysis demonstrating the shaping of functional gradients in the marmoset brain.

The discovery of functional gradients introduce a new perspective in understanding the cortical spectrum of intrinsic dynamics, as it captures major axes of functional connectivity in low-dimensional space. However, how functional gradients arise and dynamically vary remains poorly understood. In this study, we investigated the biological basis of functional gradients using awake resting-state fMRI, retrograde tracing and gene expression datasets in marmosets. We found functional gradients in marmosets showed a sensorimotor-to-visual principal gradient followed by a unimodal-to-multimodal gradient, resembling functional gradients in human children. Although strongly constrained by structural wirings, functional gradients were dynamically modulated by arousal levels. Utilizing a reduced model, we uncovered opposing effects on gradient dynamics by structural connectivity (inverted U-shape) and neuromodulatory input (U-shape) with arousal fluctuations, and dissected the contribution of individual neuromodulatory receptors. This study provides insights into biological basis of functional gradients by revealing the interaction between structural connectivity and ascending neuromodulatory system.

Aberrant brain functional and structural developments in MECP2 duplication rats.

Transgenic animal models with homologous etiology provide a promising way to pursue the neurobiological substrates of the behavioral deficits in autism spectrum disorder (ASD). Gain-of-function mutations of MECP2 cause MECP2 duplication syndrome, a severe neurological disorder with core symptoms of ASD. However, abnormal brain developments underlying the autistic-like behavioral deficits of MECP2 duplication syndrome are rarely investigated. To this end, a human MECP2 duplication (MECP2-DP) rat model was created by the bacterial artificial chromosome transgenic method. Functional and structural magnetic resonance imaging (MRI) with high-field were performed on 16 male MECP2-DP rats and 15 male wildtype rats at postnatal 28 days, 42 days, and 56 days old. Multimodal fusion analyses guided by locomotor-relevant metrics and social novelty time separately were applied to identify abnormal brain networks associated with diverse behavioral deficits induced by MECP2 duplication. Aberrant functional developments of a core network primarily composed of the dorsal medial prefrontal cortex (dmPFC) and retrosplenial cortex (RSP) were detected to associate with diverse behavioral phenotypes in MECP2-DP rats. Altered developments of gray matter volume were detected in the hippocampus and thalamus. We conclude that gain-of-function mutations of MECP2 induce aberrant functional activities in the default-mode-like network and aberrant volumetric changes in the brain, resulting in autistic-like behavioral deficits. Our results gain critical insights into the biomarker of MECP2 duplication syndrome and the neurobiological underpinnings of the behavioral deficits in ASD.

Cold-sensitive ventromedial hypothalamic neurons control homeostatic thermogenesis and social interaction-associated hyperthermia.

Homeostatic thermogenesis is an essential protective feature of endotherms. However, the specific neuronal types involved in cold-induced thermogenesis remain largely unknown. Using functional magnetic resonance imaging and in situ hybridization, we screened for cold-sensitive neurons and found preprodynorphin (PDYN)-expressing cells in the dorsal medial region of the ventromedial hypothalamus (dmVMH) to be a candidate. Subsequent in vivo calcium recording showed that cold temperature activates dmVMHPdyn neurons, whereas hot temperature suppresses them. In addition, optogenetic activation of dmVMHPdyn neurons increases the brown adipose tissue and core body temperature, heart rate, and blood pressure, whereas optogenetic inhibition shows opposite effects, supporting their role in homeostatic thermogenesis. Furthermore, we found that dmVMHPdyn neurons are linked to known thermoregulatory circuits. Importantly, dmVMHPdyn neurons also show activation during mouse social interaction, and optogenetic inhibition suppresses social interaction and associated hyperthermia. Together, our study describes dual functions of dmVMHPdyn neurons that allow coordinated regulation of body temperature and social behaviors.

Male and Female C57BL/6 Mice Respond Differently to Awake Magnetic Resonance Imaging Habituation.

Traditionally, preclinical magnetic resonance imaging (MRI) has been performed in anesthetized animals. However, anesthesia has been shown to perturb normal brain function and physiology. Such effects limit our ability to detect subtle physiological alterations in disease models and treatment studies, thus hampering discovery and compromising generality of findings. Therefore, methods for awake animal MRI are needed to study the rodent brain in its natural physiological state, free of anesthetics. Current setups for awake animal MRI rely on restraining systems to avoid animal movement during scanning. To reduce restraint stress, animals are habituated to the scanner environment prior to MRI data collection. To date, however, most awake MRI studies employ male rodents only. This is a fundamental limitation as results obtained may be pertinent only to half of the population. We characterized training and habituation responses of male and female mice to provide improved, sex-dependent training procedures for awake mouse MRI. We recorded heart rate, monitored behavioral responses (body weight and fecal boli weight) and fecal corticosterone levels (FCM) as indicators of wellbeing and stress during a 14-day progressive habituation protocol. In addition, we also assessed discomfort levels and anxiety using the mouse grimace scale (MGS) and light/dark test (LDT), respectively. All scores were compared between both groups. We found that heart rate was significantly decreased after 10 and 11 days of training for both males and females, respectively. However, the specific time course for this decrease was significantly different between males and females, and females exhibited higher anxiety levels during habituation and 14 days after habituation than males. Lastly, we also found that mean FCM levels for both groups were decreased after 11 days of MRI habituation. The present work shows that mice can be successfully trained for extended MRI sessions which is necessary for many (particularly non-fMRI) studies. Importantly, we find that males and females differ in their response to awake MRI habituation, which should be considered in future awake MRI studies that aim to include male and female mice.

Automated Skull Stripping in Mouse Functional Magnetic Resonance Imaging Analysis Using 3D U-Net.

Skull stripping is an initial and critical step in the pipeline of mouse fMRI analysis. Manual labeling of the brain usually suffers from intra- and inter-rater variability and is highly time-consuming. Hence, an automatic and efficient skull-stripping method is in high demand for mouse fMRI studies. In this study, we investigated a 3D U-Net based method for automatic brain extraction in mouse fMRI studies. Two U-Net models were separately trained on T2-weighted anatomical images and T2*-weighted functional images. The trained models were tested on both interior and exterior datasets. The 3D U-Net models yielded a higher accuracy in brain extraction from both T2-weighted images (Dice > 0.984, Jaccard index > 0.968 and Hausdorff distance < 7.7) and T2*-weighted images (Dice > 0.964, Jaccard index > 0.931 and Hausdorff distance < 3.3), compared with the two widely used mouse skull-stripping methods (RATS and SHERM). The resting-state fMRI results using automatic segmentation with the 3D U-Net models are highly consistent with those obtained by manual segmentation for both the seed-based and group independent component analysis. These results demonstrate that the 3D U-Net based method can replace manual brain extraction in mouse fMRI analysis.

Constructing the hierarchy of predictive auditory sequences in the marmoset brain.

Our brains constantly generate predictions of sensory input that are compared with actual inputs, propagate the prediction-errors through a hierarchy of brain regions, and subsequently update the internal predictions of the world. However, the essential feature of predictive coding, the notion of hierarchical depth and its neural mechanisms, remains largely unexplored. Here, we investigated the hierarchical depth of predictive auditory processing by combining functional magnetic resonance imaging (fMRI) and high-density whole-brain electrocorticography (ECoG) in marmoset monkeys during an auditory local-global paradigm in which the temporal regularities of the stimuli were designed at two hierarchical levels. The prediction-errors and prediction updates were examined as neural responses to auditory mismatches and omissions. Using fMRI, we identified a hierarchical gradient along the auditory pathway: midbrain and sensory regions represented local, shorter-time-scale predictive processing followed by associative auditory regions, whereas anterior temporal and prefrontal areas represented global, longer-time-scale sequence processing. The complementary ECoG recordings confirmed the activations at cortical surface areas and further differentiated the signals of prediction-error and update, which were transmitted via putative bottom-up γ and top-down ß oscillations, respectively. Furthermore, omission responses caused by absence of input, reflecting solely the two levels of prediction signals that are unique to the hierarchical predictive coding framework, demonstrated the hierarchical top-down process of predictions in the auditory, temporal, and prefrontal areas. Thus, our findings support the hierarchical predictive coding framework, and outline how neural networks and spatiotemporal dynamics are used to represent and arrange a hierarchical structure of auditory sequences in the marmoset brain.

Whole-brain mapping of mouse CSF flow via HEAP-METRIC phase-contrast MRI.

PURPOSE: CSF plays important roles in clearing brain waste and homeostasis. However, mapping whole-brain CSF flow in the rodents is difficult, primarily due to its assumed very low velocity. Therefore, we aimed to develop a novel phase-contrast MRI method to map whole-brain CSF flow in the mouse brain. METHODS: A novel generalized Hadamard encoding-based multi-band scheme (dubbed HEAP-METRIC, Hadamard Encoding APproach of Multi-band Excitation for short TR Imaging aCcelerating) using complex Hadamard matrix was developed and incorporated into conventional phase contrast (PC)-MRI to significantly increase SNR. RESULTS: Slow flow phantom imaging validated HEAP-METRIC PC-MRI's ability to achieve fast and accurate mapping of slow flow velocities (~102  µm/s). With the SNR gain afforded by HEAP-METRIC scheme, high-resolution (0.08 × 0.08 mm in-plane resolution and 36 0.4 mm slices) PC-MRI was completed in 21 min for whole-brain CSF flow mapping in the mouse. Using this novel method, we provide the first report of whole-brain CSF flow in the awake mouse brain with an average flow velocity of ~200 µm/s. Furthermore, HEAP-METRIC PC-MRI revealed CSF flow was reduced by isoflurane anesthesia, accompanied by reduction of glymphatic function as measured by dynamic contrast-enhanced MRI. CONCLUSION: We developed and validated a generalized HEAP-METRIC PC-MRI for mapping low velocity flow. With this method, we have achieved the first whole-brain mapping of awake mouse CSF flow and have further revealed that anesthesia reduces CSF flow velocity.

High-resolution relaxometry-based calibrated fMRI in murine brain: Metabolic differences between awake and anesthetized states.

Functional magnetic resonance imaging (fMRI) techniques using the blood-oxygen level-dependent (BOLD) signal have shown great potential as clinical biomarkers of disease. Thus, using these techniques in preclinical rodent models is an urgent need. Calibrated fMRI is a promising technique that can provide high-resolution mapping of cerebral oxygen metabolism (CMRO2). However, calibrated fMRI is difficult to use in rodent models for several reasons: rodents are anesthetized, stimulation-induced changes are small, and gas challenges induce noisy CMRO2 predictions. We used, in mice, a relaxometry-based calibrated fMRI method which uses cerebral blood flow (CBF) and the BOLD-sensitive magnetic relaxation component, R2', the same parameter derived in the deoxyhemoglobin-dilution model of calibrated fMRI. This method does not use any gas challenges, which we tested on mice in both awake and anesthetized states. As anesthesia induces a whole-brain change, our protocol allowed us to overcome the former limitations of rodent studies using calibrated fMRI. We revealed 1.5-2 times higher CMRO2, dependent upon brain region, in the awake state versus the anesthetized state. Our results agree with alternative measurements of whole-brain CMRO2 in the same mice and previous human anesthesia studies. The use of calibrated fMRI in rodents has much potential for preclinical fMRI.

Suppressing of Power Line Artifact From Electroencephalogram Measurements Using Sparsity in Frequency Domain.

Electroencephalogram (EEG) plays an important role in brain disease diagnosis and research of brain-computer interface (BCI). However, the measurements of EEG are often exposed to strong interference of power line artifact (PLA). Digital notch filters (DNFs) can be applied to remove the PLA effectively, but it also results in severe signal distortions in the time domain. To address this problem, spectrum correction (SC) based methods can be utilized. These methods estimate harmonic parameters of the PLA such that compensation signals are produced to remove the noise. In order to ensure high accuracy during harmonic parameter estimations, a novel approach is proposed in this paper. This novel approach is based on the combination of sparse representation (SR) and SC. It can deeply mine the information of PLA in the frequency domain. Firstly, a ratio-based spectrum correction (RBSC) using rectangular window is employed to make rough estimation of the harmonic parameters of PLA. Secondly, the two spectral line closest to the estimated frequency are calculated. Thirdly, the two spectral lines with high amplitudes can be utilized as input of RBSC to make finer estimations of the harmonic parameters. Finally, a compensation signal, based on the extracted harmonic parameters, is generated to suppress PLA. Numerical simulations and actual EEG signals with PLA were used to evaluate the effectiveness of the improved approach. It is verified that this approach can effectively suppress the PLA without distorting the time-domain waveform of the EEG signal.

Having Infants in the Family Group Promotes Altruistic Behavior of Marmoset Monkeys.

The common marmoset (Callithrix jacchus) has attracted much attention as a useful model for studying social behaviors [1-3]. They naturally live in a monogamous family group and exhibit cooperative breeding [4], in which parents and older siblings help to carry infants less than 2 months old [5-7]. Marmoset parents also transfer foods to their offspring, a process that may help them learn the food diet [8]. Furthermore, marmosets show spontaneous altruistic behaviors, such as providing food to non-reciprocating and genetically unrelated individuals [9]. These social habits indicate that marmosets may be a useful non-human primate model for studying parenting and altruistic behaviors, as well as underlying neural mechanisms. Using a novel rescue paradigm, we found that marmoset parents and older siblings showed strong motivation to rescue trapped young infants but not juvenile marmosets beyond 2 months of age, and infant calls alone could trigger these parents' rescue behaviors. The marmoset parents showed little rescue of each other, but young infants or infant calls could also induce such parents' mutual rescue. Moreover, all these infant- and mate-rescue behaviors depended on currently having young infants in the family group. Functional MRI studies on awake adult marmosets showed that calls from young infants, but not juvenile marmosets, elicited a large-scale activation of specific brain areas including auditory and insular cortices, and such activation was absent in marmosets not living with infants. Thus, such infant-induced modification of neural activity offers a window for examining the neural basis of altruistic behaviors in marmoset monkeys.

Correction to: Projection from the Anterior Cingulate Cortex to the Lateral Part of Mediodorsal Thalamus Modulates Vicarious Freezing Behavior.

The authors found that in this article.